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Излагается подход к расчету распространения растворенной пассивной примеси, при котором этот 
процесс рассматривается как одномерная дисперсия. Приводится способ вычисления коэффициента 
дисперсии по дисперсии скорости в створах на различных объектах, в том числе и на проточном 
Иваньковском водохранилище. 


Уравнение одномерной адвективной диспер- 
сии широко применяется при расчетах распрост- 
ранения примеси в реках и каналах [5, 10] на зна- 
чительном удалении от ее источника. Здесь рас- 
сматривается оценка коэффициента продольной 
дисперсии Д при случайном блуждании частиц 
(элементов жидкости) в потоке. 

Уравнение, описывающее распространение 
примеси. в водотоке, имеет вид [4] 
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( 1 ) 


где р - плотность вероятностей распределения 
примеси массы М вдоль оси потока, х - криволи- 
нейная ось координат (в данном случае сов- 
падающая с фарватером), V - средняя по сечению 
скорость течения. 

Решение этого уравнения в случае мгновен- 
ного точечного выпуска примеси в момент време- 
ни і = 0 в точке х 0 и при р — ► 0, если х — ► 
известно: 
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где 4*о» Г) - створ, в который попадает примесь, 
двигаясь со средней по сечению водотока скоро- 
стью, х { - координата точки наблюдения. 

Если предположить, что вся масса примеси М 
в начальный момент времени равномерно распре- 
делена в слое единичной толщины в ств. х 0 , то 
С 0 = М/5(х о), где С 0 - условная начальная концен- 
трация, 5(х) - площадь сечения водотока в ств. х. 
Учитывая, что в реальном водотоке 5(х) * сопзі, 
получаем 


С(х, г) = С 0 р(і)5(х 0 )/5( х). (3) 


Чтобы использовать уравнения (1) - (3), необ- 
ходимо соблюдать следующие условия: 


Э5 _ 2 _ ЭД 

ді дх дх 


(4) 


дѴ 

Р Тх 


~ др 


(5) 


где (2 - расход воды в сечении русла. 

В (2) Д может быть рассчитан по результатам 
измерения концентрации примеси (искусственной 
или естественной) [5], переносимой потоком, 
либо по заданным в сечении водотока распреде- 
лениям скоростей течения и коэффициентов тур- 
булентного обмена с помощью общепринятых 
формул X. Фишера и Дж. Эльдера [5, 10, 13]: 


° 0 * 0 (6) 
х [Ц(у”)-й]Ь(у п )(іу я (Іу'<1у, 
е ѵ = 0.23ЛД, 

где 11* - динамическая, ІДу) - средняя по вертикали 
скорость потока, Н(у) - глубина в точке у, Еу - коэф- 
фициент турбулентной диффузии по у. Скорость 
течения осреднена по времени и вклад турбулент- 
ной диффузии не учитывается. Чтобы воспользо- 
ваться формулами (6), необходимо задать весьма 
подробно распределения ІДу, і) и Ыу) в створе, что 
связано с техническими трудностями при натурных 
измерениях, поэтому, например, в [9] рекомендует- 
ся сочетать измерение ІДу, і) с расчетом, что позво- 
ляет уменьшить шаг интегрирования и повысить 
точность вычисления Д. 

Кроме этих двух методик при расчете Д 
часто используются эмпирические формулы [1]. 
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В большинстве случаев они могут быть представ- 
лены в виде 

О, К, В), 

где X - коэффициент гидравлического трения 
(обычно вводится через динамическую скорость 
или уклон), К - гидравлический радиус, В - шири- 
на водотока по поверхности (необязательный 
параметр). 

Набор аргументов в этой формуле вполне 
достаточен для описания равномерных потоков 
(в этом случае практически полностью определя- 
ется турбулентный обмен). В реальных реках 
дисперсия примеси может быть сильно связана с 
формой русла в сечении и в плане плохо учтенны- 
ми перечисленными параметрами. Авторы [3, 14] 
пытались учесть влияние формы русла в сечении 
и в плане, но полученные ими результаты не 
были распространены. 

Применение уравнения (1) для описания про- 
цесса дисперсии основано на аналогии между 
диффузией и дисперсией, при этом определение 
коэффициентов диффузии (или дисперсии), как 
правило, основано на градиентной теории. Как 
известно, коэффициенты диффузии могут быть 
определены с помощью теории Тейлора через 
корреляцию смещений блуждающих частиц [11]. 
Подобные соотношения могут быть получены и 
для коэффициента дисперсии. Такой подход ис- 
пользовался в [13, 14]. 

Очевидно, что частица, находящаяся при I = 0 
в точке х = 0, за время г сместится вдоль траекто- 
рии по оси х на расстояние 

г 
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Л = -Р-, X = [ Ц(х)сІХ , х=\ й(х)сІх' , (7) 
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заменим в предыдущей формуле переменные 
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Осредним обе части полученного соотношения 
по сечению в ств. х: 
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Здесь ( ) означают осреднение по сечению 
потока: 
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где ІГ = і/ - Ц, 5 = х - х, г(%) - аналог автокорре- 
ляционной функции, <?„(*) и о и (х) - среднеквадра- 
тические отклонения скорости воды от 0(х) и 
0{х) в сечениях іих соответственно. 


Домножая обе части равенства на [/(/) и учи- 
тывая, что 

...... 1 (Іх 2 
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1 Положив мгновенные значения скорости и расстояния 
равными і) + и и х + х‘ соответственно, после осреднения 
по вероятности получим 

о 


+ 1 н'(г)м'(т)4т. 
о 

Здесь первый интеграл определяет адвекцию примеси в 
потоке, а второй - турбулентную диффузию. 


Пусть 0(х) ~ 0(х) и с и (х ) » о и (х), т.е. скорость 
изменяется незначительно при значениях І*-*! 
приблизительно до 2 Тогда 


1 А (х + а 2 ) 

2 (И 


<*«(*) 

ад 



т, 

о 


где а] - дисперсия распределения примеси вдоль 
оси х. 


Так как 
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и, если при х — *- «> интеграл от г(Е,) в (10) стре- 
мится к конечному пределу Ь, 


Я. 



(И) 


Такое представление накладывает определенные 

ограничения на поле скорости: при 7/ — ► 0 
выражение (11) теряет смысл. 

Параметр ^ является масштабом дисперсии и 
может быть истолкован как минимальное рассто- 
яние между створами, начиная с которого распре- 
деления скоростей в этих створах можно считать 
независимыми. Вычислить г(^) непосредственно 
невозможно, поэтому естественно попытаться 
связать величину 7, с одним из масштабов рус- 
ловых форм, например со средней протяженнос- 
тью излучин или радиусом корреляции автокор- 
реляционной функции отклонений тальвега от 
условий оси [2]. В общем случае масштаб I і может 
зависеть от многих параметров (например, от 
фрактальной размерности [7]), в том числе харак- 
теризующих картину русла в плане. 

В [2] измерялись отклонения фарватера Ни- 
жней Волги от условной оси и вычислялась их 
автокорреляционная функция г(^) вдоль этой оси. 
Соответствующий интегральный масштаб, рас- 
считанный до первого нуля г(^), оказался прибли- 
женно равным 5 В (В - средняя ширина русла по 
поверхности). Средняя протяженность излучин 
равнинных рек обычно оценивается [7] величи- 
ной (5 - 7) В. Поскольку поперечный профиль 
скорости на протяжении излучины полностью пе- 
рестраивается, выбор ее длины в качестве мас- 
штаба дисперсии кажется оправданным. При 
формальном переходе к прямому руслу в этом 
случае 7, — ► а г(^) — ► сопзі, так как распреде- 
ление скорости во всех сечениях одинаково. Для 
принятых условий из (10) для прямого канала по- 
лучим 




т 


(11а) 


При ламинарном течении г - 1. Физическая 
определенность такого перехода также служит 
доводом в пользу выбора масштаба Ь, равного 
средней длине излучины. 

В дальнейшем выражение 7,-5 В принято в 
качестве ориентировочного для равнинных рек с 
ограниченным меандрированием. Однако такой 
способ вычисления Ь необходимо проверить на 
экспериментальном материале. 

Дисперсию скорости в створе, также входя- 
щую в (11), можно определить по формуле 

В Ку) 

= Ш (і2) 

0 0 


ИЛИ 

Вк(у) 

о]= Д Г [ Щх, у) - ІДу)] г Ат.<1у + 

0 ° в (12а) 

+ Діуоо-йѴ 

о 

В (12а) первый интеграл определяет диспер- 
сию распределения скорости по вертикали, вто- 
рой — по ширине. 

Для прямой экспериментальной проверки при- 
веденной методики необходимо располагать ре- 
зультатами измерений поля скорости по сечени- 
ям, коэффициента дисперсии, полученного с по- 
мощью трассеров, и достаточно подробным 
планом ряда рек и каналов. К сожалению, такими 
данными мы не располагаем, поэтому контроль- 
ные расчеты для рассматриваемых объектов осу- 
ществлялись с помощью (6), если это было воз- 
можно, и эмпирических формул [1]. 

По данным, приведенным в [1], оценим измен- 
чивость отношения ЦВ, предположив, что рас- 
сматриваемые потоки плоские с профилем ско- 
рости 

ІЛу,і) = (1 + а)ІЛу)[і/И(у)]« (13) 

где а - показатель степени, который соответству- 
ет суммарному в сечении реки трению и опреде- 
ляется зависимостью [8] 

а = 7/. /к 7/, (14) 

к = 0.4 - постоянная Кармана. 

С помощью соотношений (12) - (14) получим 

сі = (а 2 &Ѵ(1+2а), Ь = О х /а г и . (15) 

Таким образом, 7,/в ([1], табл. 2) изменяется 
приблизительно от 4 до 20, причем для прямых лот- 
ков -от 17 до 20. В этой таблице Д т /7(7/* е [9.8; 558]. 

Далее был выбран объект (р. Педедзе [6]) с по- 
дробно промеренными в сечении полем скорости 
и профилем дна, для которого 7) г рассчитывался 
по (6) при 

б у = 0.23/?7/.. (16) 

С помощью измеренных значений ІДу, і) методом 
наименьших квадратов определялся показатель 
степени а (в логарифмическом приближении) в 
предположении справедливости закона распреде- 
ления скорости (13). Динамическая скорость мо- 
жет быть определена по формуле (14), в которой 
использовалось среднеквадратическое значение а 
(осреднение осуществлялось по ширине реки). 
В результате было получено, что О х = 13.3 м 2 /с. 
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С, г/м 3 



Изменение во времени концентрации фосфора в 
ств. д. Плоски. Точки - эксперимент, кривая - расчет. 


С помощью этих данных вычисления осуще- 
ствлялись также по эмпирической формуле 


КУ. 


0.88 



(в) 


и ' 


(17) 


рекомендованной в [1] при В - 60 - 70 м. В резуль- 
тате й х - 13.1 м 2 /с. 


циент 


Вычисленный по (1 1), (12) при Ь = 5В коэффи- 
;нт дисперсии 2 в этом случае равен 12 м 2 /с. 


Как видно, результаты трех вариантов расчета 
О г весьма близки. 


По результатам измерений, проведенных на 
р. Киржач (подробное описание эксперимента 
дано в [10, гл. 5]), О х рассчитывался по (1 1), (12) и 
( 1 7) при I, = 5 В. Полученные значения коэффици- 
ентов дисперсии 2.9 и 3.5 м 2 /с с учетом ограничен- 
ности применения эмпирических формул и услов- 
ности выбора Ь удовлетворительно согласуются 
одно с другим. Не исключено также, что в этом 
случае Ь > 5 В. 

Для Иваньковского водохранилища расчеты 
й х осуществлялись по результатам измерения 
расходов в восьми створах: г. Дубна ( х = 0), 
о. Клинцы (л: = 10.2 км), г. Корчева' ( х = 18.4 км), 
г. Конаково (. х = 36.6 км), д. Плоски (х = 42.4 км), 
о. Низовка ( х = 53.7 км), с. Городня (.г = 70.4 км), 
г. Эммаус (л: = 93.4 км). 

Расстояния определены от ств. Дубна по фар- 
ватеру с помощью лоцманской карты. 


На участке г. Эммаус-с. Городня практически 
реализуется режим речного течения, поэтому 
здесь можно воспользоваться эмпирической фор- 
мулой (для широких русел) [1] 


О х 

ки. 


= 2138 


А 

; 


1.63 


(1В) 


2 При вычислении интеграла (6) использовалась сплайн-ин- 
терполяция профилей скорости и дна. Шаг интегрирова- 
ния, обеспечивающий необходимую точность, <В/ 200, т.е. 
существенно меньше рекомендованного в [8] значения 
В/ 20. Интеграл (12) вычислялся методом трапеций. 

3 Город полностью затоплен. 


При определении II* по таблицам выбирался ко- 
эффициент шероховатости (0.03 - 0.06), а рас- 
чет осуществлялся по формулам Маннинга е 
Ш ези [12]. Если принять, что среднее для указан- 
ного интервала значение коэффициента шерохо- 
ватости равно 0.045, то й х = 7.9 м 2 /с. 

Согласно материалам Конаковского краевед- 
ческого музея, ширина русла реки в межень в кон- 
це XIX в. на участке выше с. Городня составляла 
-160 м, т.е. Ь - 800 м при Ь/В = 5. В расчетах эта 
величина была оценочной, так как естественно 
предположить, что русло и пойма реки сформи- 
ровались до создания водохранилища. 

Коэффициент дисперсии в ств. г. Эммаус, 
вычисленный по (11), (12) при Ь = 800 м, равен 
8.5 м 2 /с (среднее значение). 

Аналогично (по (1 1), (12) при Ь = 800 м) значе- 
ние й х было вычислено и в других створах. Раз- 
брос результатов этих вычислений оказался зна- 
чительным. Очевидно, это связано с возникнове- 
нием в водохранилище циркуляционных течений, 
ускоряющих перемешивание водных масс и уве- 
личивающих дисперсию с\ Наибольший разброс 
Ц т получился в створах о. Клинцы и г. Корчева 
(наиболее сложный широкий участок водохрани- 
лища). В ств. г. Дубна разброс несколько уменьша- 
ется (сказывается организующее влияние водо- 
слива). Очевидно, что наиболее простая картина 
течений реализуется при проточном установив- 
шемся режиме течения, чему должно соответст- 
вовать минимальное значение дисперсии скоро- 
сти и, следовательно, Г>,. 

Коэффициенты дисперсии, отобранные в со- 
ответствии с этим принципом, имеют следующие 
значения: 10.6 (г. Дубна), 11.1 (о. Клинцы), 
17 (г. Корчева), 18.0 (г. Конаково), 15.5 (д. Плоски), 
19 (о. Низовка), 12 м 2 /с (с. Городня). Среднее по 
длине водохранилище (г. Дубна-г. Эммаус) значе- 
ние О х равно 14 м 2 /с. 

Еще одна оценка Ц т была получена по наблю- 
давшемуся в апреле в 1984 г. в ств. д. Плоски крат- 
ковременному всплеску концентрации общего 
фосфора 4 (рисунок). Предполагалось, что этот 
импульс сформировался в Шошинском заливе в 
результате начавшегося паводка на реках Шоше 
и Ламе, текущих с юга, и его дисперсия намного 
меньше дисперсии распределения, наблюдавше- 
гося в ств. д. Плоски. 

Оценка была выполнена по (2), (3) методом 
наименьших квадратов. При этом величина 
С(х, I) - С тіп нормировалась на максимальное на- 
блюдавшееся значение (д:, = .ѵ(х 0 г)) и результат 
логарифмировался (С^ = 0.021 г/м 3 - фоновая 
концентрация фосфора). 

4 Учитывая значительное время добегания и низкую темпе- 
ратуру воды (~ 0°С), при вычислении й х можно прибли- 
женно считать фосфор нереагирующей примесью. 
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Влияние распределения скорости течения II 
вдоль водохранилища учитывалось следующим 
образом. 

При установившемся режиме и в отсутствие 
сильного ветра поле скоростей Иваньковского 
водохранилища определяется расходом воды 
через плотину (). В силу малой боковой приточно- 
сти (<10%) на рассматриваемом участке этот рас- 
ход можно считать независимым от х, в то же вре- 
мя скорость течения II плавно убывает при при- 
ближении к плотине. 

По данным измерений расходов в названных 
выше створах была подобрана функция 


&(х, О) = а 0 + а 1 х + а 2 () + а рс 2 + а л & + а (19) 

где х - расстояние по фарватеру, отсчитываемое 
от ств. г. Дубна. 

Если 0(х, О) выражено в см/с, х - в км, а (2 - в 
м 3 /с, то а 0 = 2.423, а х = -3.790 х ІО" 2 ; а 2 = -6.630 х 
х ІО" 3 , а ъ = 9.296 х 10Л а 4 = 2.469 х 10" 5 , а 5 = 1.978 х 
х КГ 4 . 

Аппроксимация (19) удовлетворительна (сред- 
неквадратическое отклонение для всех створов 
составляет 1.35 см/с) при () е [60, 800]. Так как 
(2(х) = 0(х, О) 5(х, 0, то выражение (19) опреде- 
ляет также и сглаженное значение 5(х, 0, необ- 
ходимое для вычисления концентрации по (3): 


5(х 0 )/5(х) = 0 (х)/ 0 (х 0 ). 


Среднее время добегания между створами х 0 и 
х, определяется соотношением 


или с учетом (19) 


к 


СІХ 


Щх, 0 


2 2а 3 х + а х + а 5 <2 
1 = — агсІ§ т= 

-/А Л 


( 20 ) 


где А = 4 а 3 (а 0 + а 2 () + а^О 2 ) - {а х + а 5 0 2 . 

Примесь, выпущенная при Г = 0 в ств. х 0 , дости- 
гает за время I створа 


, Ді§(? - г, Д/2) - (а, + а 5 0 

Фо >0 = 2аз ’ (21) 

где і - агсі§[(2азХ 0 + а х + а 5 { 2)1 ^ ]. 

Коэффициент дисперсии (по распределению 
концентрации фосфора), вычисленный с по- 
мощью рассмотренной методики с применением 
(2), (3), (21), равен 13.6 м 2 /с. Условный центр мас- 
сы растворенного фосфора при I = 0 считался 
расположенным в 65 км от ств. г. Дубна, а <2 = 


= 150 м 3 /с. Кривая, рассчитанная по (2), (3), (21), 
приведена на рисунке. 

Оценка выполнения условий (5) при средних 

по х значениях производных, р и II при удалении 
от источника на 25 км показала, что правая часть 
неравенства в несколько раз превышает левую 
(при больших удалениях от источника это отно- 
шение увеличивается). Если для всего водохра- 
нилища среднее значение Ѳ х равно 14 м 2 /с, что 
хорошо согласуется с расчетом коэффициента 
дисперсии по распределению концентрации рас- 
творенного фосфора, то при рассматриваемом 
режиме течения выполняются и условия (4). 

Таким образом, метод определения коэффи- 
циента продольной дисперсии по дисперсии ско- 
рости в створах дает в рассмотренных случаях 
удовлетворительный результат. 

Авторы благодарят В.Д. Казмирука за по- 
мощь в работе над статьей. 
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